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Abstract 

Strain redistribution corrections were devel- 
oped for a simplified inelastic an’ysis procedure 
to economically calculate material cyclic response 
at the critical location of a structure for life 
prediction purposes. The method was based on the 
assumption that the plastic region in the structure 
is local and the total strain history required for 
input can be defined from elastic finite-element 
analyses. Cyclic stress-strain behavior was repre- 
sented by a bil near kinematic hardening model. The 
simplified procedure has been found to predict 
stress-strain response with reasonable accuracy for 
thermally cycled problems but needs improvement for 
mechanically load-cycled problems. 

This study derived and incorporated Neuber-type 
corrections in the simplified procedure to account 
for local total strain redistribution under cyclic 
mechanical loading. The corrected simplified method 
was exercised on a mechanically load-cycled bench- 
mark notched plate problem. Excellent agreement 
was found between the predicted material response 
and nonlinear finite-element solutions for the 
problem. The simplified analysis computer program 
used 0.3 percent of the CPU time required for a 
nonlinear finite-element analysis. 

Introduction 

Nonlinear, finite-element computer codes are 
generally too costly to use as a design tool for 
hot section components of aircraft gas turbine 
engines. Computing costs are further Increased by 
the presence of high thermal gradients and geomet- 
rical irregularities, such as cooling holes. Most 
hot section components are of such geometrical 
complexity as to necessitate three-dimensional 
analyses, frequently with substructuring. Three- 
dimensional, nonlinear finite-element analyses are 
prohibitively time consuming and expensive to con- 
duct in the early design stages for combustor and 
turbine structures. To improve the durability of 
these components, it is necessary to develop 
simpler and more economical methods for represent- 
ing the structural response of materials under 
cyclic loading. 


♦Similar material will be presented in Local 
Strain Redistribution Corrections for a Simplified 
Inelastic Analysis Procedure Based on an Elastic 
Finite - Analysis by Albert Kaufman and Shol Y. 
Hwang (NASA TP-2421). 


A program has been underway at NASA Lewis to 
develop a simplified procedure 1- '’ for performing 
nonlinear structural analysis using only an elastic 
finite element solution or strain gage data as 
input. Development of the simplified method was 
based on the assumption that the inelastic regions 
in the structure are constrained by the surrounding 
elastic material. This implies that the total 
strain history can be defined by an elastic analy- 
sis. Initial development of the method di^not 
account for any strain redistribution. A computer 
program (ANSYMP) was created to predict the stress- 
strain history at the critical fatigue location of 
a thermomechanical ly cycled structure from elastic 
input data. Appropriate material stress-strain and 
creep properties and plasticity hardening models 
were incorporated into the program. Effective 
stresses and equivalent plastic strains are approx- 
imated by an iterative and incremental solution 
procedure. Creep effects can be obtained on the 
basis of stress relaxation at constant strain, 
cumulative creep at constant stress or a combina- 
tion of stress relaxation and creep accumulation. 

The simplified procedure has been found to pre- 
dict critical location stress-strain response with 
reasonable accuracy relative to nonlinear finite- 
element analyses for thermally cycled problems. 
However, the limitation of no strain redistribution 
is most likely to be violated in the case of 
mechanically loaded structures, especially in the 
vicinity of stress concentrations. Nonlinear 
finite-element analyses of notched plate specimens 
subjected to cyclic mechanical loading have shown 
that the total strain range at the critical loca- 
tion can be significantly larger than would be 
predicted from elastic solutions. 

This Sfudy derived and incorporated corrections 
in the simplified procedure to account for local 
total strain redistribution and residual stresses 
due to mechanical load cycling. These corrections 
would remove some of the limitations and extend the 
applicability of the procedure. The corrections 
were based on the Neuber rule 6 relating the theo- 
retical stress concentration factor to the actual 
stress and strain concentration ^actors. 

Two variations of a benchmark notched plate 
problem 6 were analytically examined to verify the 
accuracy of the improved simplified method. Veri- 
fication was made through comparison with three- 
dimensional nonlinear, finite-element analyses. 
Cyclic stress-strain and creep properties for 
Inconel 718 alloy, a kinematic hardening model and 
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the von Mises yield criterion were used for the 
benchmark notch specimen problems. Elastic and 
elastic-plastic finite element analyses were per- 
formed using the MARC nonlinear finite-element 
computer code.' The elastic solutions for the 
critical locations were used as input data for the 
simplified analysis computer code. Verification 
of the improved simplified procedure was made on 
the basis of how well it was able to duplicate the 
stress-strain hysteresis loops from MARC elastic- 
plastic analyses of the benchmark notch problem. 

Nomenclature 

e nominal total strain at net section 

E modulus of elasticity 

E e modified modulus of elasticity 

Ep work hardening slope (Fig. 1) 

K,n temperature-dependent constants in cyclic 

stress-strain equation 

K t theoretical elastic stress or strain con- 
centration factor 

K c actual strain concentration factor 

K 0 actual stress concentration factor 

S effective nominal elastic stress at nth 

increment 

S* effective local elastic stress (real or 
imaginary) 

Sj maximum effective local elastic stress 

c equivalent total strain, unmodified 

e* equivalent total strain, modified 

e e uncorrected equivalent elastic strain 

cp equivalent plastic strain 

cp maximum equivalent plastic strain in cycle 

(Fig. 1) 

ep equivalent local residual strain 

a effective stress (loading), unmodified 

o' maximum effective stress in cycle (Fig. 1) 

02 effective stress (unloading), unmodified 

o* effective stress (loading), modified 

o| effective stress (unloading), modified 

o m maximum effective stress, unmodified 

om maximum effective stress, modified 

or local residual stress 

oy} initial yield stress in loading part of 

cycle 


oyi initial yield stress in unloading part of 

cycle 

u Poisson's coefficient 

Analytical Procedure 

Improvements were made to a simplified ine- 
lastic procedure for calculating the stress-strain 
history at the critical fatigue location of a 
structure subjected to cyclic thermomechanical 
loading. The basic assumption ir, the initial 
development of the procedure was that the inelastic 
region is local and constrained from redistribution 
by the surrounding elastic material. It follows 
from this assumption that the total strain history 
at the critical location can be defined by an 
elastic solution. Justification for the assumption 
of elastic constraint of loc^l inelasticity can be 
found in finite element analyses of thermally 
cycled structures, such as comoustor liners, 
aircooled turbine blades, and wedge fatigue speci- 
mens which have shown that the total strain ranges 
from elastic and nonlinear solutions are in close 
agreement. However, this assumption is invalid for 
structures subjected mainly to mechanical load 
cycling, as was shown* for the benchmark notch 
problem. 

In this study, a strain redistribution factor 
(SRF) was derived to account for local total strain 
redistribution under applied cyclic loads. This 
correction also accounted for residual strains 
induced by plastic yielding. The strain redistri- 
bution factor is applied to the ideal local total 
strain obtained from the elastic solution. These 
corrections were incorporated in a versio - of the 
ANSYMP program using a kinematic hardening model 
to characterize the yield surface under cycling. 

A representation of a cyclic stress-strain curve 
by a bilinear kinematic hardening model is illus- 
trated in Fig. 1. The loci of the tips of the 
cyclic curves is described by the relation 

O = K( c p) n (1) 

The work hardening slope, Ep, for the kinematic 
hardening model was determined from energy consid- 
erations to give the same strain energy, as indi- 
cated by the enclosed area in Fig. 1, as the actual 
stress-strain curve. This work hardening slope 
will be defined by 



The von Mises criterion was used, thereby con- 
verting the total strain from a uniaxial stress- 
strain curve into a modified equivalent total 
strain. I’ this modified elastic part of the 
equivalent total strain corresponds to the measured 
elastic strain multiplied by 2(1 ♦ u)/3. This 
relationship must be taken into account in a multi- 
axial stress-strain fielo when applying results 
from elastic finite-element analyses. As a result, 
all stresses and strains from the elastic solution 
are expressed in terms of von Mises effective 
stress and modified equivalent total strain; these 
values arc assigned signs based on the Dominant 
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principal stress and strain directions. For com- 
putational convenience, the factor 2(1 ♦ u ) / 3 was 
Included in the elastic modulus and the modified 
elastic modulus will be designated by E e where 
E e - 3/ ( 2( 1 ♦ w) )E. 

The elastic input data are subdivided into a 
sufficient number of increments to define the 
stress-strain cycle. Dwell times are specified for 
increments which require creep analysis. The 
increments are analyzed sequentially to obtain the 
cumulative plastic and creep strains and to track 
the yield surface. An iterative process is used 
to calculate the yield stresses for increments 
undergoing plastic straining. First, an estimated 
plastic strain is assumed for calculating an ini- 
tial yield stress from the stress-strain properties 
and the simulated hardening model. Second, a new 
plastic strain is calculated as the difference 
between the totai strain and the elastic and creep 
strain components. The yield stress is then 
recalculated using the new plastic strain. This 
Iterative process is repeated until the new and 
previous plastic strains agree within a tolerance 
of 1 percent. 

Neuber's rule for approximating local stresses 
and strains in the plastic region was used as the 
basis for the development of the strain redistri- 
bution corrections. With minor modifications, this 
rule provides acceptable predictions of local 
stresses and strains in the loading part of a 
stable cycle and less accurate predictions in the 
unloading part of the cycle. The relative failure 
of Neuber's rule for the unloading part of the 
cycle has been mainly due to neglect of the local 
residual stresses. Regardless of the part of the 
cycle being studied, an adjustment must be made to 
the stresses to account for the work hardening 
slope of the cyclic stress-plastic strain curve. 

The strain redistribution analysis for the 
stable cycle is divided into three stages; elastic- 
plastic loading, elastic unloading and plastic 
unloading. Figure 2 is presented to further define 
the symbols used in the governing equations for 
each of these stages. 

Elastic-Plastic Loadi ng 


0 “ (* ~ j ® y1 * E p c when c > ^ (6) 

Combining Eqs. (4) to (6) gives the local 
strain as 



(7) 


where 


A . 2E P 



The local stresses and strains calculat'd by 
using the work hardening slope in Eqs. (6) and (7) 
have to be modified to account for the elastic 
strain change between o and oyi. The modified 
local stress, o* , corresponding to the local 
strain, t*, can be expressed as 



and the modified local strain, c*, as 


C 


* 


C 



(9) 


The plastic strain, c p at the stress level o* 
is 


( 10 ) 


To relate the ideal local inelastic strain to 
the local total strain, a strain redistribution 
factor (SRF) is introduced. The SRF at the nth 
increment on the stress-strain curve can be writ- 
ten as 


Neuber's rule specifies a relationship between 
the theoretical stress concentration factor, K^, 

.ind the actual stress concentration factor, K 0 , 
ind the actual strain concentration factor, K t , 
of the form 

K. 2 . K K (3) 

t o c 

Substituting local and net section stresses and, 
strains into Eq. (3) gives 



where Kt also equals S*/S or o « S*^/(cE e ). 

The cyclic stress-strain relation is assumed to 
take the form 

o » E g c when c < o^/ E g (5) 


( SRF ) n “ *n ” r~ <“> 

e 

To improve the accuracy of the predicted local 
stress-strain history, a strain redistribution 
factor, SRF, was generated for each increment dur- 
ing the loading part of tha cycle. This SFR 
represents the cumulative total strain correction 
up to the nth point on the stress-strain curve. 
Therefore, the total strain at that increment is 
the uncorrected elastic strain, c e , for that 
increment plus the sum of all the incremental SRFs 
up to that increment. 

n 

'n-- e „*Z[ (SRF) t-l SRF >M] <“) 

1-1 

Equation (12) is another way of expressing c* 
given in Eq. (9). 


or 
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Elastic Unloading 

The local stress at the nth increment in the 
elastic unloading part of the cycle (Fig. 2) is 
simply 


o n - o* - (S* - S*) 
n m m n 


and the corresponding local total strain is 

* . 
rn 

T 


E_ ■ E* - 

n m 


S* - S* 
m n 


(13) 


(14) 


B - 4E n 


■ p \ , + S * 2 

C/°yl e R £7 ' °R e R 


The actual local stresses and strains are 
obtained by a series of iterations between the 
modified and unmodified values. Two iterations, 
as shown below, have been found sufficient to give 
an acceptably convergent solution. Iterations do 
not appear to be necessary for the loading part of 
the cycle. The tpodified local stress, 02 , cor- 
responding to t is expressed as 


[(SRF). - (SRF) 1 _J (15) 


or expressed In a form similar to Eq. (12) 
n 
c 

o en 

i-1 

For an isothermal problem where an incremental 
approach is not required, £■;. (15) can be simpli- 
fied to t n *■ cgn + ( SRF ) m . 

Inelastic Unloading 

Upon complete remova. of the applied loaa, 
residual stresses remain across the net section. 
Denoting the local residual stress as or, then 


- s: 

m m 


(16) 


and the complimentary local residual strain, Er, is 

(17) 


R 


Again the cyclic stress strain relation in the 
plastic region is assumed to take a linear form 




+ E P c 


(18) 


As a kinematic hardening model is used in the 
analysis, the new initial yield stress is given by 


yi 


a* — 2a 
m yi 


(19) 


Considering the zero stress and lastic strain 
at o m the new origin, as shown in Fig. 2, the 
local stress and total strain can be expressed as 


J 2 m \r; °yi 


♦ i - 




where 03 is obtained from the following steps 
1 * £ p‘ 




S| ■ t — 


(° ~ 0 *) 

En 


°1 = l 1 " ET ) °yi + E P C 1 




C ? = E - 

E„ 


(Oj- O^) 

' "E'n 


J 2 "l 1 '^i°yi T t P c 2 


°yi + E P £ 2 

and the modified loral total strain, e*, as 

C* = t 


0 “ °2 °2 - °2 . 
r— *p 


(23) 

> 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 


K S + 


( 20 ) 


The plastic strain, r p , at the stress level 
uj is 


and 


e = K e ♦ c 


( 21 ) 


Combining Eqs. (T), (18), (20), and (21) leads 
to the solution for 


. .-A- V(A 2 ♦ B) 


( 22 ) 


where 


v°yi 


'P*R 


£p = — [ (32) 

e 

The strain redistribution factor at the nth 
increment on the stress-strain curve is 

/ mr \ °n _ °n2 0 n2 -0 n2 S *r\ + °R 

(SRF) n ’ * n r n 1 P (33) 


T 


t: 


or, in terms of the uncorrected elastic strain, the 
total strain at the nth increment can be written as 


c* * r ♦ 
n en 


(SRF) m ‘ X [ (SRF) i ' (SRF, i-ll (34) 

i-1 
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where 


(SRF) m 



Equation (34) Is simply another form of expres- 
sing r* previously given In Eq. (31). 

Another version of the ANSYMP computer program 
was created to implement the improved simplified 
analytical procedure. The program consists of the 
main executive routine, ANSYMP, and four subrou- 
tines, ELAS, YIELD, CREEP, and SHIFT. The incre- 
mental elastic data and temperatures are read into 
subroutine ELAS. Material stress-strain properties 
as a function of temperature and a simulated hard- 
ening model are incorporated in subroutine YIELD 
and the creep characteristics are incorporated in 
subroutine CREEP. Subroutine SHIFT is required to 
update the temperature effects on the yield stress 
shift. SHIFT also serves the function of deciding 
the future direction of the yield surface under 
nonisothermal conditions by determining the rela- 
tion of future to past themal loading. 

The calculational scheme initially follows the 
effective stress-equivalent strain input data from 
subroutine ELAS until the occurrence of initial 
yielding. The stress-strain solution then proceeds 
along the yield surface as determined from the 
stress-strain properties in subroutine YIELD At 
each increment during yielding the stress shift 
(difference between new yield stress and stress 
predicted from elastic analysis) from the original 
input data is calculated. Elastic load reversal 
is signaled when the input stress is less than the 
yield stress from the previous increment. During 
elastic unloading, the stresses are translated from 
the original elastic analysis solution by the 
amount of the calculated stress shift. Reverse 
yielding occurs when the stress reaches the reverse 
yield surface as determined from the hardening 
model incorporated in subroutine YIELD. Again, the 
solution follows the yield surface until another 
I ad reversal is indicated when the stress based 
on the shifted elastic solution is less than the 
yield stress. The elastic response during load 
reversal is obtained by translating the original 
elastic solution according to the new stress shift 
calculated during reversed yielding. The stress- 
strain response for subsequent cycles is computed 
by repeating this procedure of identifying load 
reversals, tracking reverse yield surfaces and 
translating the original elastic solution during 
elastic loading and unloading. Creep computations 
are performed for increments involving dwell times 
using the creep equation and strain hardening rule 
inccrporated in subroutine CREEP. Depending on the 
natjre of the problem, the creep effects are deter- 
mined on the basis of one of the three options pro- 
vided in the suoroutine. 

The code will automatically avoid the Neuber- 
type strain redistribution corrections for thermal 
loading problems where there are no applied mechan- 
ical loads. Without applied loads, the Neuber 
method would be inapplicable since Eq. (4) would 
have zero net stresses and strains in the denomi- 
nator. Provision is also made for the user to 
circumvent the corrections for other situations 
where they should not be taken into account. These 
situations include locally strain controlled prob- 
lems and, as will be shown, problems where the 
total strain input is based on strain measurements 
rather than elastic finite-element analyses. 


Since the stable cyclic stress-strain curve is 
a function oi the plastic strain range, it Is nec- 
essary to iterate between the maximum plastic 
strain used to specify the stress-strain curve and 
the analytical solution. In both the MARC code and 
the original ANSYMP code, this required rerunning 
the problem until the assumed and calculated plas- 
tic strains were in reasonable agreement. In the 
improved version of the ANSYMP program, the Itera- 
tion was automated so that the analysis was 
repeated using the previously calculated maximum 
plastic strain to define the stress-strain curve. 
This iterative process is continued until the 
specified and calculated maximum plastic strains 
agree within 10 percent, usually within three 
iterations. 

The computer program was verified by conducting 
simplified analyses for the oenchmark notch problem 
and comparing the results to those from MARC non- 
linear analyses. The geometry of the benchmark 
notch specimen is illustrated in Fig. 3. This 
specimen was tested under Isothermal conditions as 
part of a program to provide controlled strain data 
for constitutive model verification.” A MARC 
analysis of this problem using kinematic hardening 
demonstrated excellent agreement with experimental 
data. 4 Two variations of this problem were ana- 
lyzed; one with the initial loading carried out to 
a plastic strain of 0.4 percent and the second with 
the specimen loaded to a maximum plastic strain of 
0.6 percent with a small amount of strain ratchet- 
ting. A kinematic hardening model was used with 
cyclic stress-strain data for Inconel 718 alloy.® 
Nonlinear and elastic MARC analyses of this problem 
were performed using approximately 600 triangular 
plane strain elements to model a quarter segment 
bounded by planes of symmetry as shown in Fig. 4. 

The MARC solutions shown for the benchmark notch 
specimen were computed at the closest Gaussian 
integration point to the root of the notch. 

Discussion of Analytical Results 

The results of the improved simplified inelas- 
tic analyses of the benchmarx notch problems is 
discussed herein. Comparisons are made with MARC 
inelastic solutions. Stress-strain cycles used for 
comparison purposes ars in terms of effective 
stresses and equivalent total strains using the 
von Mises yield cr<*crion computed with signs based 
on the assigned signs in the original elastic solu- 
tion. The entire discussion is based on the crit- 
ical location at ti.o notch root. 

The benchmark notch test was conducted by 
mechanical load cycling at a constant temperature 
of 649 *C. A mechanically load’d structure, espe- 
cially wherz the peak strain occurs at a stress 
raiser, is most likely to violate the basic assump- 
tion of the simplified approach that strain redis- 
tribution is prevented by containment of the local 
plastic region by the surrounding elastic material. 
It was shown! that the total strain range from 
the MARC elastic-plastic analysis was 20 percent 
greater than that obtained from the MARC elastic 
analysis. This foreshortening of the elastic 
strain range caused the simplified procedure to 
truncate the stress-strain hysteresis loop as shown 
in Fig. 5(a). When the input total strain history 
was based on optical strain measurements at the 
notch root, the agreement between the simplified 
and MARC elastic-plastic stress-strain hysteresis 
loops was good as demonstrated in Fig. 5(b). Both 
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the ANSYMP and MARC elastic-plastic analyses gave 
stable stress-strain hysteresis loops for the sec- 
ond cycle. It Is to be noted that when the simpli- 
fied analysis uses strain measurements rather than 
elastic finite-element solutions, strain redistri- 
bution corrections are unnecessary and should not 
be applied. 

In Fig. 5(c) similar comparisons are shown 
using the simplified procedure corrected for strain 
redistribution on the loading part of the cycle and 
the MARC elastic solution as Input. Neuber correc- 
tions for residual stresses and strain redistribu- 
tion on the unloading part of the cycle were not 
Implemented for this case. The truncation of the 
stress-strain hysteresis loop shown In Fig. 5(a) 
was eliminated and the A 1YMP solution showed good 
agreement with the MARC solution. This demon- 
strates the significant Improvement In accuracy 
that can be attained by applying the Neuber correc- 
tion In the plastic region, even without the resid- 
ual stress correction during unloading. 

When both the loading and unloading strain 
redistribution corrections were applied, the agree- 
ment between the predicted and MARC elastic-plastic 
results was even better, as shown In Fig. 5(d). 

The original benchmark problem analyzed in 
Fig. 5, was a completely closed cycle and the max- 
imum plastic strain was about 0.4 percent. To 
exercise the improved simplified procedure on an 
even more severe case, the mechanical loading was 
increased so that the plastic strain reached 
approximately 0.6 percent and strain ratchetting 
was Induced. Again the agreement between the 
ANSYMP and MARC inelastic solutions was excellent, 
as shown in Fig. 6. It should be noted that the 
simplified procedure was able to capture the strain 
ratchetting on the second cycle. 

The slight discrepancies between the the ANSYMP 
and MARC hysteresis loops in Figs. 5(d) and 6 were 
due primarily to the iteration process built into 
ANSYMP which resulted in the use of a cyclic 
stress-strain curve slightly different from that 
used In the MARC elastic-plastic analysis. The 
ANSYMP analyses of the benchmark notch problem used 

0.3 percent of the central processor unit (CPU) 
time required by the MARC nonlinear analyses; this 
Is without considering that the latter had to be 
run several times to obtain a cyclic stress-strain 
curve compatible with the calculated maximum 
plastic strain. 

Summary of Results 

An improved simplified analysis procedure was 
developed for calculating the stress-strain history 
at the critical location of a thermomechanical ly 
cycled structure. This Improved procedure incor- 
porated Neuber-type corrections to account for 
strain redistribution and residual stresses due to 
plastic strain reversals. Analytical predictions 
from the simplified procedure for benchmark notched 
plate problems were compared to nonlinear finite- 
element solutions. The following general conclu- 
sions were drawn from the evaluation of the method: 

1. The predicted stress-strain response based 
on elastic finite-element solutions with the 
Neuber-type corrections showed excellent agreement 
with elastic-plastic finite-element solutions using 
the MARC program. 


2. The predicted stress-strain response using 
only the corrections for plastic yielding on the 
loading part of the cycle showed good agreement 
with the MARC elastic-plastic solution. The cor- 
rections for the unloading part of the cycle were 
of secondary importaice. 

3. The strain redistribution corrections should 
ot be applied when the input total strain history 

is based on local strain measurements rather than 
elastic finite-element analyses. It is also known 
from previous evaluations of the simplified proce- 
dure that these corrections should not be applied 
in thermal loading problems. For cases where a 
structure is subjected to a combination of thermal 
and mechanical loads, a method will have to be 
developed to partition the thermal and mechanical 
stresses so that the strain redistribution correc- 
tions are only applied to the mechanical stresses. 

4. Stress-strain hysteresis loops were computed 
at the critical location oi the structure using 

0.3 percent of the CPU time required for elastic- 
plastic finite-element analyses. 
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(b) Loading part of eyrie. 
Figure 2. - Description of symbols. 
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(c) Unloading part of cycle. 
Figure 2. - Concluded. 


71.12 (RADIUS)- 


11 (RADIUSH \ 


L- 

1 1 ' v 

2.54 

__ 



_J ✓ 



t ^ 1 

i r 1 

__ 



4.902 


17.958- 
24. 13 


/< 


i i 
i i 


SYMMETRIC 
CENTERLINE 
- 19.05 
(RADIUS) 


i 



I 



[ ur\ — i 

t 

1 


254 00 


19. 05 (DIAM) 


25. 4WDIAM) 


Figure 3. - Benchmark notch specimen (kj * 1.9). (All dimen- 
sions in mm.) 
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Figure 4 - Benchmark specimen finite-element model. 
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(b) Strain measurements as input, uncoriected for 
strain redistribution. 

Figure 5. - Comparison of ANSYMP and MARC results 
for benchmark problem (0. 4 * maximum plastic strain). 
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Figure 6. - Comparison of ANSYMP and MARC results for 
benchmark problem (0. 6% maximum plastic strain). 


